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Abstract: We amalgamate the many experimental limits on the abb coupling of a light 
CP-odd Higgs boson, a, including model-dependence coming from the ratio of the ati to the 
abb coupling. We then employ these limits to analyze the extent to which a light a can make 
a significant contribution to the discrepancy, Aa^, between the experimentally observed 
a^ and that predicted by the standard model. In a "model-independent" framework and 
in the context of a general two-Higgs-doublet model this is a significant possibility. In 
contrast, the minimal supersymmetric model is too strongly constrained (after combining 
experimental and theoretical input) to allow a CP-odd-a explanation of Aa^. The next- 
to-minimal supersymmetric model allows more freedom and the light a of the model could 
explain the full Aa^ if 9.2 GeV < rua < 12 GeV, or contribute substantially for larger m^, 
if tan f3 is large. 



There have been numerous studies fl], |, |, |, |, |, 0, ||, |, |To|, [n|, |l|] of the extent to 
which the Higgs sector could contribute to the anomalous magnetic moment of the muon, 
a^, with current focus on whether it could be used to explain some portion of the now ~ 3a 
positive deviation of with respect to the Standard Model (SM) prediction. The numerical 
deviation is variously quoted as Aa^ ~ (27.5 ±8.4) x 10-^° or (27.7 ±9.3) x 10~i° |2|. 
It is becoming increasingly likely that this deviation can only be explained by new physics 
of some kind and a beyond-the-standard-model Higgs sector has always been a prime 
candidate. 

Precision electroweak data and direct LEP limits on a light CP-even scalar suggest that 
it should have SM-like couplings and substantial mass, in which case its contribution to 
will only be of order few x 10^^^. Thus, we will focus on the possible contribution, 6a^, of 
a light CP-odd Higgs boson, a, of a CP-conserving Higgs sector, for which it is critical |l| 



to include the two- loop Barr-Zee type diagrams |14| since the one- loop contribution is 
negative whereas the two-loop contribution is positive in popular models. 

Of particular interest is the rua < 2mh region, for which a light Higgs, h, with SM-like 
WW, ZZ and fermionic couplings can have mass ruh ~ 100 GeV while still being consistent 
with LEP data by virtue of /i —> aa —> 4r decays being dominant |15, |l6|, (see also 

1 19, ^). Such a Higgs provides perfect agreement with the rather compelling precision 
electroweak constraints, and for BR[h aa) ^ 0.75 also provides an explanation for the 
~ 2.3(7 excess observed at LEP in e^e~ Zbb in the region M^-^ ~ 100 GeV. We term 
this the "ideal" Higgs scenario. More generally, we will only consider models for which the 
ZZh coupling is SM-like (implying zero Zha coupling and therefore no lower limits on rua 
coming from e^e^ — s- ha at LEP) and nih is such as to give good agreement with precision 
electroweak data. 

Possible contributions to by the a depend crucially on the afi^j^i^, abb and ati 
couplings defined via 

We assume a Higgs model in which Ca^-^+ = Car-T+ = C'abfe' typified by a two- 
Higgs-doublet model (2HDM) of either type-I or type-II (a 2HDM contains Higgs bosons 
h, H with mu > TUh, a and /i^), or more generally if the lepton and quark masses are 
generated by the same combination of Higgs fields. (Much larger values of relative to 
those we find below are possible in models in which r = {Caf^-^+ = CaT-T+)/C^fjj^ ^» 1 — 
such models include those in which the muon and tau masses are generated by different 
Higgs fields than the b mass. For r ^ 1, our results for should be rescaled by r.) 
In a 2HDM of type-II and in the MSSM, 0^^-^+ = C^^^ = tan/3 (where tan/3 = hu/h^ 
is the ratio of the vacuum expectation values for the doublets giving mass to up-type 
quarks vs. down-type quarks) and C^^j = cot (3. In the NMSSM the expressions for 
Cafj.-fi+ = C'afefc ^^'^ ^att include an additional factor discussed later. In a type-I 2HDM, 
Cat,-ti+ = Car-T+ = C'afeb = ~^att = " cot /3. In the most general Higgs model, (7^^-^+, 
Car-T+ ) C'^feb ^att ^^^^ more complicated functions of the vevs of the Higgs fields and 
the structure of the Yukawa couplings. In this paper, we assume Cq^-^+ = (7^^-^+ = C^^^ 
but allow for general values of R'^^^ = C^f^j^/C^^j. We consider only positive values of R'^j^ 
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since only these are of relevance for explaining the observed positive Aa^ and positive 
values are typical of most models. 
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Figure 1: Upper limit, C'^^, on \C^g^^ as a function of ma coming directly from experimental 
data. In the case of limits based on a — ^ t+t^, curves for i?^/^ — 0.5 (red), 1 (blue), 2 (green), 10 
(black) are shown with i?b/t = 10 giving the lowest curves and i?^/^ = 0.5 giving the highest curves. 



Limits on \C^^^\ come from T — s- 7a decays at B factories and e^e~ bba production 
at LEP. For nia < 2mj- the strongest limits come from the old (90% CL - all other limits 

'yX), where X is assumed to 
limits on T — > 7a ^ "yrr are 



employed here are 95% CL) CUSB-II limits ^ on BR{T 
be visible. For 2mr < rua < 9.2 GeV, the recent CLEO-III ^ 

the strongest (in interpreting these limits one must account for the value of BR{a t ' t 
— this in turn depends on Rh/t, but very weakly for i?^,/^ > 2, see below). For 9.2 GeV < 
iria < mx, mixing of the a with various rjh and xo bound states becomes crucial p3|| . 
Ref. gives results for C™|^ in this nia range without taking this mixing into account 
but notes that their limits cannot be relied upon for nia > 9-2 GeV. Whether additional 
limits can be extracted from lepton non-universality studies in the 9.2 < rUa < ttt-t region is 
being studied [^]. OPAL limits |^5[ (which assume BR{a — > r+r^) = 1) on e^e~ bbrr 
become numerically relevant for roughly 9 GeV < nia < 2nif,. Ref. [|25| converts these 



limits to limits on the abb coupling using the modeling of [23|. These are the only limits in 
the nij- < nia < 2ni{y range and continue to be relevant up to 12 GeV. Above nia = 2nif, 
these abb coupling limits become quite weak due to the rjh — a mixing and the decrease 
of BR{a — > rr). For nia ^ 12 GeV, limits on the abb coupling can be extracted from 
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e^e^ bba bbbb |2^. The maximum value of \C^fjj^\ allowed by all these various limits, 
C™|^, is shown in Fig. |l| as a function of rria for several values of Rf,/t {Rb/t = 0-5, 1, 2, 10). 
Note that there is almost no dependence of C"^^ on R^^n for i?;,/^ > 2. Values of |C^j,^| 
above 50 raise issues of non-perturbativity of the abb coupling and are likely to be in 
conflict with Tevatron limits on bba production [^]. C"^^ depends on ii;,/^ when the 
CLEO-III T — 5- 7a — > 7T^r~ or OPAL bba — s- bbr^r^ limits are the most relevant. What 
is new in this paper is the systematic incorporation of the R^/t dependence of and 
the systematic incorporation of the C^^^ limits in the context of predictions for 5a^ in a 
wide class of models. 

In the case of the simple 2HDM(II), where C^^j = Rf,/t = tan/3, values of rua for 
which tan/3 > C™v^(tan/3) are not allowed in the model context. These disallowed regions 

abb ^ ' ' " 

typically emerge in the range rUa < 8 GeV for tan/? = 1 rising to rUa ^ 10 GeV for 
higher tan/?; at higher tan/3 values they have a complicated structure that we will discuss 
later. In addition, a disallowed region also arises over a limited nia range starting from 
rUa > 12 GeV when tan/3 > 18, the larger the value of tan/3 the larger the interval. For 
example, for tan/3 = 50 the DELPHI limits imply that the 2HDM(II) is not consistent for 
12 < rUa < 37 GeV and the OPAL and Upsilon limits imply that the 2HDM(II) is not 
consistent for rUa < 10 GeV. In contrast, for tan/3 = 10 the 2HDM(II) model is always 
consistent with the DELPHI limits and is only inconsistent (with CLEO-III and CUSB 
limits) for nia ^ 9 GeV. These 2HDM(II) results are an update of the results obtained in 
Q. The results in all other models, in particular in the NMSSM context are new. ^ 

We will now explore the implications for a^. Since the two-loop contributions include 
that with a t-loop as well as those with b and r loops, we must specify the value of C^^j 
relative to C^j,^ in order to compute the contribution of a to for a given C^j^^ value. In a 
2HDM of type-II, including the MSSM and NMSSM, C^^j = cot /3 and after including the 
two-loop diagrams (5a^ > for ma > 2.6,2,0 GeV if tan/3 > 5,3,1. In a type-I 2HDM, 
^aii — ~^abb ~ '^^^ f^' Then, the (dominant) top-loop Barr-Zee type diagram gives a 
negative contribution to a^ and 5a ^ is negative for all m^. Only models with positive R^n 
are of relevance for explaining the observed positive Aa^. Results for 6a ^ employing the 
^abb ~ ^Tbl^ limits as a function of nia and and taking Ri,/i = 1, 3, 10 and 50 are 
plotted in Fig. (For iJf,/^ < 1, simply multiply the Rh/t = 1 curve by l/R^^^.) 

To a good approximation, > 50 is equivalent to dropping the two- loop diagram 
containing the top quark and gives the smallest result. Since (for positive C'^b^/C'^tj) the 
two-loop top diagram enters with the same (positive) sign as the b and r two-loop diagrams, 
the largest 5a^ values are obtained for the smallest R^/t when using upper limits on the abb 
coupling as input. As a result, we see in Fig. |2| that for lower Ri,/t values (1 < R^/t ^ 3) 
any value of rua ^ 9 GeV would make it possible to obtain 6a^ = IS.a^ ~ 27.5 x 10"^'' 
for some choice of C^^^ < C™^. For ii^/^ < 0.2, for which C^^j enters non-perturbative 
territory, 6a^ > Aa^ if C^^^ = for all nia so that agreement could always be obtained 

for some C^^-j^ < C^^^. However, for i?^/^ > 10 the full discrepancy can only be explained 
if 10 GeV < < 12 GeV or nia ^ 36 GeV. Recall, however, that the value of C^f" 

^Several months after arXiv submission of this paper, similar results for the NMSSM were obtained in 

ii- 
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Figure 2: The value of 5a^ from CP-odd a loops is plotted as a function of ttIq for i?^/^ ~ 1 (black, 
solid), = 3 (green, dashes), R^/t = 10 (blue, dots) and i?^/^ = 50 (red, long dash, short dot, 
lowest curve), assuming maximal abb coupling, C™|'^, from Fig. |^. Also shown as the dotdash cyan 
curve is the 2HDM(II) prediction for tan/3 = C™|'^(tan/3), i.e. the largest possible (self-consistent) 
choice of tan/3 within the 2HDM(II) context. For rua ^ 9 GeV this latter curve coincides with the 
Rh/t = 50 curve. 

extracted from the data in the former region relies on the modeling for the a — rjb mixing 
employed in the experimental analysis. Also, for nia > 36 GeV and R(,/f > 10, (5a^ = Aa^ 
requires non-perturbative C^^^ > 50 values. 

Of course, it is interesting to know what value of C^g^ < C^^^ is needed in order to 
match the observed Aa^ = 27.5 x 10"^'^ for those and Rh/t values for which this is 
possible. The results for the general case in which C^f^j^ is not correlated with fif,/^ are 
plotted in Fig. ^ In general, for low values of Rjj/i (for which the top loop is a major 
contributor to 5a ^) rather modest values of C^^^ will reproduce the observed Aa^. As Rb/t 
increases, the bottom loop diagrams must reproduce Aa^ on their own and increasingly 
large values of C^^j are required. As we shall see, one particularly interesting range of ttt-q 
for Rf)/f > 10 is 9.9 GeV < nia ^ 12 GeV. In Fig. ^, we see that in this range the 
observed Aa^, = 27.5 x 10-^° is matched for C^^^ in the range 28 < C^^^ < 32 for R^/t > 10 
when 9.9 GeV < nia < 12 GeV. 

The above results are modified in the context of more restrictive models. Fig. |^ shows 
the results for 5a^ in the type-II 2HDM, in which C^^,^ = Ri^ji = tan /?, obtained for various 
tan (5 values. In the type-II 2HDM, the value of 5a^ is determined once tan (3 and ma 
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Figure 3: The value of of C^j^j required in order that (Ja^ ~ 27.5 x 10^^" is plotted as a function 
of nia for _R(,/t = 1 (black, solid), i?;,/f =3 (green, dashes), — 10 (blue, dots) and i?;,/^ = 50 
(red, long dash, short dot, highest curve), for those choices of ma such that the required C^f,^ is less 
than as plotted in Fig. |]. Gaps for any given Rhjt curve correspond to TOq values for which 

^abb ^^^^ would be required. 



are specified. Unlike the very general case just considered, for which i?^/^ is not related to 
C^j^j, in the 2HDM(II) context one cannot have large C^i^ without having large R^jt^ which 
then minimizes the very important (positive) top loop contribution. Thus, the largest (5a^ 
values are now obtained with large tan/? values. The possibilities are also constrained by 
the requirement that tan/3 cannot exceed C™5^(tan/3). The gaps in the curves of Fig. ^ 
are those regions where tan/3 > C™i^^(tan /?). The result is that in order to obtain a value 

abb ^ ' 

of ba^ of order 27.5 x 10^^'^ that also has tan/3 < C™|^(tan/3) requires a rather precisely 
fixed value of tan/3 ~ 30 — 32 and rua ~ 9.9 — 12 GeV (see the tan/3 = 32 dotdash cyan 
curve). In the context of the most general CP-conserving type-II 2HDM, any value in the 
above small range is not excluded using combined Zh and ha LEP data |29| so long as 
vrih ^ 60 GeV; and, there are no limits on if vrih ^ 100 GeV. Further, contributions to 
the precision electroweak observables S and T are tiny if mn = 171^+ when h has SM-like 
ZZh coupling. As a further remark, we note from trends as tan /? increases apparent in 
Fig. ^ (lower plot) that for tan (3 values above 50 (i.e. outside the perturbative limit on 
this coupling) one will not be able to have tan/3 < C™|^(tan/?) in the rUa < 12 GeV zone, 
but that at some largish value of nia above about 40 GeV one will be able to achieve a 
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Figure 4: The value of from CP-odd a loops is plotted as a function of iria for tan /? = 1 
(black, solid), tan/3 = 3 (green, dashes), tan/? — 10 (blue, dots), tan/3 = 22 (black, dotdash), 
tan/3 = 32 (cyan, dotdash) and tan/3 = 50 (red, long dash, short dot, highest curve), assuming 
the 2HDM(II) model with i?f,/j — tan/3 and requiring that tan/3 < C™|'^(tan/3). Omitted regions 
are those for which tan/3 > C™|''(tan/3) as plotted in Fig. |l|. Note the multiple gaps for the 
tan/3 = 22, 32, 50 cases in the 10 GeV < nia < 11 GeV region. An intersection of the solid red line 
with a 2HDM(II) curve essentially only occurs in the tan (3 — 32 case. 



at (5a^ = 27.5 x 10 



match to Aa^. This is because the DELPHI limits on C^^j^ deteriorate so rapidly as rua 
increases above 40 GeV. 
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As a further perspective on the 2HDM(II) results, we plot in Fig. |2| the largest possible 
value of 5a^ within the 2HDM(II) as a function of rua (the dotdash cyan curve). This 
maximal value is obtained when tan/3 = C™5^(tan/3) (i.e. for the largest self-consistent 
choice of tan P such that C^^,^ = tan /3). Again it is apparent that 5a^ can match (or exceed) 
27.5 X 10~^'' in the range 9.9 ^ rUa ^ 12 GeV. And, to repeat, matching in this range is 
always achieved for tan/3 ~ 30 — 32. 

The ability to achieve (5a^ = Aa^ is much more constrained in the popular Minimal 
Supersymmetric Model (MSSM). In the MSSM, the LEP lower limit on nia is of order 



90 — 100 GeV, depending upon tan/3 and precise model inputs |^. For rUa > 90 GeV, 
dUfj, = Ao^ is only achievable for C^^^j = tan /3 well above the upper bound of 50 employed 
here. (Of course, if the MSSM sparticles are light, their contributions could yield the 
observed Aa^ ||3l[| .) 

The Next-to-Minimal Supersymmetric model (NMSSM) provides more fertile ground. 
The NMSSM is obtained by adding a singlet superfield S to the MSSM Higgs superfields 
Hu and Hd- Ref. |32] was the first to consider the NMSSM Higgs sector phenomenology in 



detail. The scalar component of S contains one CP-even and one CP-odd scalar field. The 
resulting Higgs sector thus contains three CP-even Higgs bosons (/ii,2,3) and two CP-odd 
Higgs bosons (01,2)) all of which can have a singlet component. A convenient program 
for exploring the NMSSM Higgs sector is NMHDECAY [||, ||]. We will not consider 
contributions to from sparticles as recently studied in [0, 12 1. 



The NMSSM is especially attractive in that it allows for the "ideal" Higgs sector 
described earlier with 171^^ ~ 100 GeV, consistent with LEP data if nia^ < 2m;, and 
BR{hi aici) > 0.75. For nia^ > 2mfe, one must have m/i^ > 110 GeV to avoid LEP 
bounds. (But, for 110 GeV ^ mh^ ^ 163 GeV, so long as the ZZh\ coupling is SM-like 
the agreement with precision electroweak data is still within the 95% CL limit unless only 
the "leptonic" determination of sin^ is employed in the precision electroweak analysis; 
the latter yields a much higher CL for the overall fit and requires m/^j < 105 GeV at 95% 
CL — see for details). 

The most crucial parameter for the NMSSM analysis is cos Qa defined by 

fli = cos Oaumssm + sin 9Aas, (2) 

where umssm is the CP-odd (doublet) scalar in the MSSM sector of the NMSSM and as is 
the additional CP-odd singlet scalar of the NMSSM. In terms of cos^^^, Cq^^—^+ — ^abb — 
cos 9a tan /3 and C^fi = cos 9a cot (3. 

Before proceeding, we consider possible constraints from precision electroweak data. 
Since the light SM-like hi already gives good agreement, the rest of the Higgs sector should 
give a small contribution to S and T (assuming sparticle contributions are not substantial). 
One finds that if rua^ is in the range considered and hi is SM-like, then it is typically the 
case that either /12 or /13 is mainly singlet, denoted /ig, and the other, denoted here as ho, 
is mainly doublet. Further, the Zhgai coupling is very tiny while the Zh^ai coupling is 
maximal and mi^+ ~ ~ m/i^. With these inputs, one finds that the extra contributions 
from the Higgs sector to S and T are very small and the excellent agreement with precision 
electroweak constraints coming from the hi is preserved. 
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Figure 5: cos 6'™^'^ in the NMSSM (where — cos 6*^ tan/3) as a function of vfia after requiring 
cos 6'™"^ tan /? = using the C'™|'^ values of Fig. |l| The different curves correspond to tan l3 — 1 

(black, soHd, upper curve), 3 (green, dashes), 10 (blue, dots), 32 (cyan, dotdash) and 50 (red, long 
dash, short dot, lowest curve). 



Let us now consider oi = a contributions to a^j for various fixed tan /3 values. Then, 
cos 6a is constrained by the requirement that C^^^ = cos 6a tan (3 < C^^^^, which constrains 
cos 6a to very small values for low and large tan (3. However, no matter what the value 
of tan P, the extra freedom of adjusting cos 6a does allow us to avoid gaps in for which 
^abb ^ ^^l^' '^^i^' ™ turn, will give us more possibilities for (5a^. Inputting the values of 
^max ^ function of TUa we obtain the results of Fig. H for the maximum allowed value of 

abb '-' 

COS 6 A as a function of nia for various tan/3 values. 

We now turn to the resulting NMSSM predictions for a^. The value of 5a^ is largest 
for cos 6a = cos 6^^^^. The resulting values of <5a^ are plotted as a function of nia in Fig. |6|. 
As in the generic case, the strong constraints from Upsilon physics imply that significant 
contributions to are not possible until nia exceeds roughly 9.2 GeV. To understand why 
6a increases with increasing tan/3 for nia > 12 GeV, whereas it decreases with increasing 
tan/3 for low nia, we first note that the 2-loop, top-loop contribution to 5a^ is independent 
of tan/? (because of a Ca^- ^+Catt structure that is tan /3- independent), whereas the 2-loop 
bottom-loop contribution increases as tan^ (3 (because of a Ca^- ^+C^^-^ oc tan^ /3 structure). 
Numerically, before including the extra tan^ /3 factor for the latter, the 2-loop, top-loop 
contribution is much larger than the 2-loop, bottom-loop contribution. Of course, both 
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Figure 6: The maximum contribution of the CP-odd a to as a function of iria in the NMSSM, for 
which C^jj; — cos 9 A tan/3, after employing cos 6 a = cosd'^'^^ where cos 9'^'^^ is plotted for different 
tan /3 values in Fig. ||. Curve notation is as in Fig. ||. The horizontal solid red line is located at 
da^ = 27.5 X 10-^°. 



contributions are multiplied by (cos^a)^- Thus, when is independent of tan/? and 

cos^™''^ = 1 (as for nia > 12 GeV and tan/3 < 20) the resulting 5a^ will always increase 
with tan/3. However, at low nia, the very strong Upsilon constraints on C^^^jj imply that 
cos^™^^ rapidly decreases with increasing tan/3 which suppresses the numerically more 
important 2-loop, top-loop contribution resulting in smaller (5a^ as tan/3 increases. 

From Fig. ^, we observe that the maximal 6a^ can exceed Aa^ = 27.5 x lO^^*' for 
9.9 GeV ^ rua ^ 12 GeV if tan /3 > 32, with an almost precise match to this value of Aa^ for 
tan /3 = 32 (or for tan /3 as low as tan /? = 30 — see the 2HDM discussion). For tan /3 = 50, 
one can match Aa^ by using a value of cos 0a below cos 6'™'^^. (As discussed below, the 
fact that matching is possible for 9.9 GeV ^ rUa ^ 2ms is particularly interesting in the 
context of the ideal Higgs scenario.) Further, the maximal 6a^ is in the 7 — 20 x lO"^'' 
range for 12 GeV < < 48 GeV for tan/3 = 32 and for 12 GeV < < 70 GeV for 
tan/3 = 50. 

At this point, it is worth stressing the other desirable features of the m/j ~ 100 GeV, 
nia ^ 2mB, BR{h aa) > 0.75 scenario as discussed in |]l5|, |l^, 18 1. These references 
examined the degree to which obtaining the observed value of mz requires very precisely 
tuned values of the GUT scale parameters of the MSSM and NMSSM. One finds that in 
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any supersymmetric model this finetuning is always minimized for GUT scale parameters 
that yield a SM-like h with nih < 100 GeV, something that is only consistent with LEP 
data if the h has unexpected decays that reduce the h ^ bb branching ratio while not 
contributing to h ^ bbbb (also strongly constrained by LEP data). A Higgs sector with a 
light a for which BR{h — > aa) > 0.75 and with small enough that a decays to BB final 
states are disallowed (i.e. < 10.56 GeV) provides a very natural possibility for allowing 
minimal finetuning. The NMSSM provides one possible example. 

In conclusion, the combined limits from T decays and bba Yukawa production at LEP, 
along with the perturbativity requirement of C^^^ < 50, imply that the entire discrep- 
ancy of Aa^ 30 X 10^^'^ cannot have a purely Higgs sector explanation without going 
beyond the MSSM. In the less-constrained NMSSM, achieving 5a^ ~ Aa^ requires rela- 
tively high tan (5 and a value of ma between about 10 GeV and 2mB- On the one hand, this 
is a highly motivated ma region in the NMSSM since, as described earlier, it would allow an 
"ideal" SM-like h with mh ^ 100 GeV decaying mainly via h ^ aa —> 4r. Such an h would 
escape LEP limits while allowing for low m^-finetuning. However, on the other hand, in 
the NMSSM ma < 2mB most naturally arises when close to the U{1)r symmetry limit. In 
this case, the a is mainly singlet, implying that cos 9 a is small and that C^^^ = cos 9 a tan /3 
is typically 0{1) |17], whereas C^^^ ~ 30 is needed to match the observed Aa^. 



Nonetheless, the possibility that a CP-odd a with 10 GeV < ma < 12 GeV could ex- 
plain the anomaly should be taken seriously. Thus, finding techniques to experimentally 
probe for an a in the 10 GeV < < 12 GeV region should be a high priority. Such new 
techniques could either end up limiting C^^^ sufficiently that Aa^ cannot be explained in 
the 2HDM(II) or NMSSM frameworks or else actually allow a discovery of a light a. Of 
course, this is a region in which tji, — a mixing will surely be a complication. 

As an aside, one must not forget that in supersymmetric models sparticle loops could 
have two important roles: (i) they could directly yield large contributions to a^; and (ii) 
they could modify the relations between Ca^-^+-, Car-r+i ^ahb ^"^^ ^att- 

If one goes beyond the MSSM and NMSSM Higgs sectors to the more general type-II 
2HDM, then, keeping C^^u < 50, only an a with 10 GeV < < 12 GeV with C^^^ ~ 30-32 
could give 6a^ = Aa^. (A type-I 2HDM gives negative (5a^ that is large for > 8 GeV 
^abb ~ ^Sf^ therefore strongly disfavored by the observed positive Aa^.) 

Obtaining the observed Aa^ in the most general Higgs model for which the abb coupling 
magnitude is disconnected from the ratio i?^^^ of the abb to ati couplings is generically 
possible so long as R'^^^ > 0. For i?^/^ = 1, ma > 8 GeV and a relatively modest value 
of Cjjj,^ (well below the maximum allowed) will yield 5a^ = Aa^. As -R^/j increases, the 
required C^^j; increases. For larger Rh/t, there are regions of ma for which the required C^j,^ 
exceeds the upper experimental bound, C^|^- Further, 5a^ = Aa^ cannot be achieved 
above an i^h/^-dependent maximum if C^f,^ < 50 is imposed. For Rjj/f < 0.2, even very 
low values of ma will yield the observed Aa^ for an appropriate choice of C^^^ < C™^. 
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